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SUMMARY 


A simplified lifting-surface theory is applied to the problem of 
evaluating span loading due to flap deflection for arbitrary wing plan 
forms. With the resulting procedure, the effects of flap deflection an 
the span loading and associated aerodynamic' characteristics can be easily 
computed for any wing which is symmetrical about the root chord and which 
has a straight quarter-chord line over the wing semispan. The effects of 
compressibility and spanwise variation of section lift— curve slope are 
taken into account by the procedure. 

For the case of straight— tapered wings, the load distribution and 
lift due to flap deflection are presented in chart form for a broad range 
of wing plan farms. Experimental values of flap effectiveness are compared 
with the values obtained from these charts for several wings. The compari- 
son shows good agreement between theory and experiment. 

/ 

The method presented can also be used to calculate the downwash in the 
vertical center of the wake of a wing which has arbitrary spanwise loading. 


INTRODUCTION 


Reference 1 presents a procedure which is based on simplified lifting- 
surface theory and which permits the rapid determination of the spanwise 
distribution of symmetric lift on wings having essentially arbitrary plan 
form. With reference 1, the symmetric spanwise loading can be found far 
any symmetric angle— of-attack distribution such as that due to spanwise 
continuous twist or camber . The method satisfies the boundary condition 
at. only a finite number of spanwise points; thus, to insure accuracy for 
the loading obtained by this method, the angle— of-attack distribution 
must be continuous and must not vary excessively. If the angle— of-attack 
distribution is discontinuous, such as that due to partial wing— 6pan flaps, 
or varies excessively, the boundary condition should be satisfied > at suffi- 
cient points to insure a correct answer. This normally entails an excessive 





2 


NA.CA TN 2278 


amount of labor, since the computations required increase approximately 
as the square of the number of points used. 

An alternative procedure is to use a method which satisfies the 
boundary condition at only a few points where this boundary condition is 
adjusted so that at these points the calculated loading is identical with 
that which would be predicted by a method satisfying the boundary condi- 
tion at all points. This specific boundary condition for the case of 
antisymmetrically deflected ailerons is developed in reference 2, which 
is based on the zero-aspect— ratio theory of reference 3. A similar devel- 
opment can determine the specific boundary condition for the case of sym- 
metrically deflected flaps. 

It is the purpose of the present report to determine the specific 
boundary condition and to provide simple methods for computing the aero- 
dynamic characteristics due to flap deflection. 

' NOTATION 


Vn 


aspect ratio 


/V\ 

V s / 


coefficient depending on wing geometry and indicating the influence 
of symmetric loading at span station n on the downwash angle at 
span station v 


span of the wing measured perpendicular to the plane of symmetry, 
feet 


av 


c l 

C T 


wing chord, 1 feet 

mean wing chord 1 (jO > feet 


local lift coefficient 


f local lift 


lift coefficient 


V <l c 

f total lift 


V 


qS 


rate of change of lift coefficient with flap deflection, per radian 


c L c av 


rate of change of lift coefficient with wing angle of attack, per 
radian 

/ \ 

2AG N 


spanwise loading coefficient for unit lift coefficient ( J 

^ °L / 


■‘•Measured parallel to the plane of symmetry. 




m arbitrary number of span stations 

q free-stream dynamic pressure, pounds per square foot 

5 wing area, square feet 

V free— stream velocity, feet per second 

w induced velocity, normal to the lifting surface, positive for down- 
wash, feet per second 

x longitudinal coordinate measured from the lateral plane through 
the quarter chord of the wing— root chord, 1 feet 

y lateral coordinate measured from the wing root perpendicular to the 
plane of symmetry, feet 

a wing angle of attack, 1 radians 

a v section angle of attack at span station v, 1 radians 

— lift-effectiveness parameter 1 
dS 

0 compressibility parameter 
F circulation, feet squared per second 

6 angle of deflection of flap, 1 radians 
1 See footnote 1, page 2. 
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angle of deflection of flap measured perpendicular to the hinge 
line, radians 

dimensionless lateral coordinate 

dimensionless flap span on one wing panel, measured perpendicular 
to the plane of symmetry, from the wing root outboard for inhoard 
flaps , and from the wing tip inboard for outboard flaps 



trigonometric spanwise coordinate defining the spanwise end of the 
flap, radians 

ratio of section lift-curve slope at span station v to — , 
both at the same Mach number 1 ^ 



A 

X 

I 

<P 


sweep angle of the wing quarter— chord line, positive for sweepback, 
degrees 

sweep parameter which includes the effects of compressibility- 
^tan — 1 , degrees 

wing taper ratio ( tip 

\root chord/ 

dimensionless longitudinal coordinate 


(4) 


b/: 

trigonometric spanwise coordinate (cos -1 tj), radians 


Subscripts 


n 


k 

f 


integer pertaining to span station associated with the loading 
rut 


G ■ 003 Et) ' 


integer pertaining to span station associated with the downwash 




= cos 


V7T 

m+1 


half— integer pertaining to span station ( T| = cos 
pertaining to flaps 


kit 

m+1 J 


See footnote 1, page 2. 
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t wing tip 

r wing root 

av average or mean 

1 denoting full-wing chord flaps ( — = 1 

V c ‘ , 

DEVELOPMENT OF METHOD 


The simplified lifting-surface theory used herein replaces a lifting 
surface hy a lifting vortex located at the wing quarter-chord line. The 
boundary condition specified for determining the vortex strength distri- 
bution is that, along the three-quarter-chord line of the wing, there shall 
be no flow through the lifting surface. 

The mathematical method of solution and discussion of the limitatio n s 
of the theory are given in detail in MCA Rep. 921, 1948 (reference 1). 

The same mathematical procedures and limitations apply to the present 
report. 

The mathematical solution consists of satisfying the boundary condi- 
tion at a finite number of points spanwise along the three— quarter— chord 
line. Then, if the boundary condition (i.e., the spanwise -twist along 
the three— quarter— chord line) varies smoothly and is continuous, relatively 
few points that satisfy the condition are needed and computations can be . 
completed without excessive labor. Thus it has been fo und that only 
seven spanwise points which satisfy the boundary condition are needed for 
accurate computations of the loading due to symmetric continuous twist 
(reference 1) and that due to antisymmetric continuous twist (reference 2). 

For the case of irregular and discontinuous distribution of twist 
such as that required to simulate the span loading due to deflection of 
a partial— span flap, it becomes necessary to satisfy the boundary condi- 
tion at a large number of points, particularly near the point of discon- 
tinuity, in order that stiff icient accuracy will be insured. .Since the 
computation labor increases approximately as the square of the number of 
points considered, the foregoing method becomes less practical. An alter- 
nate approach is to alter the boundary condition for the method involving 
few points at which the boundary condition is satisfied, such that at 
these points the loading is identical with that obtained by the method 
for which the boundary condition is satisfied at all points. This process 
is practical because the spanwise loading of wings becomes independent of 
plan form as aspect ratio approaches zero (reference 3)» Thus the solu- 
tions for the zero-aspect— ratio wing which satisfy the boundary condition 
at all points can be used, to obtain the proper adjustment of the boundary 
condition for a method that satisfies the boundary" condition at only a 
few points. With the altered boundary condition, the method satisfying 
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conditions at only a few points will give accurate results for plan 
forms of very low aspect ratio, and should give good accuracy for wings 
of high aspect ratio since the alteration of the boundary condition can 
be kept small by proper choice of flap spans, and also since loading 
distribution is affected less by plan -form than by the boundary condi- 
tion. The mathematical procedure for determining the altered boundary 
conditions is developed in appendix A. 

Although an altered boundary condition can be found for an arbitrary 
position of the point of discontinuity, the alteration can be made a 
minimum and of small order by locating the discontinuity at the mean 
position between two points at which the boundary condition is satisfied. 
With the point of discontinuity defined, the flap span becomes defined. 
The results for arbitrary flap spans are found by cross plots of results 
for the defined spans. 

When the flap spans and corresponding twists (altered boundary con- 
ditions) have been determined, the method is applied in the same manner 
as that of reference 1. In this way, (m+l)/2 linear equations in terms 
of symmetric— loading distribution are obtained which satisfy the wing 
angle— of-attack conditions 2 at the three— quarter— chord line at m 
stations, where m is ah arbitrary odd integer. The application in the 
present report is with m = 7* The linear equations are represented by 
the summations 


n=i 

where the integer subscripts v and n pertain to* span stations 
rj = cos (Vrt/8) and (cos nit/8) , respectively. To obtain the loading 
coefficients G n from the simultaneous solution of equation (l), it 
remains only to evaluate the coefficients ay n (which depend on wing 
geometry, compressibility, and section lift-curve slope) and the spanwise 
variation of the symmetric angle of, attack cty. The reader should note 
that with a v constant (spanwise) equation (l) represents additional 
loading, and that with ay varying symmetrically and continuously it 
represents basic loading. 

Evaluation of a., 

■ vn 

As shown in reference 1, a Vn can be conviently presented as a 

function of two parameters, 3 namely, a sweep parameter defined as 
Ap = tan -1 [(tan A) /p] and a wing chord distribution parameter Hy 1 
defined by 

2 As discussed in reference 1, aT is not limited to small values but can 
be as large an angle as desired, provided separation does not occur. 
3 The effects of compressibility and section lift-curve slope are equiva- 
lent to a change in wing plan forms and can be taken into account by a 
proper adjustment of the a vn values. 
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( 2 ) 


where 

k v ratio of experimental section lift— curve slope at span station v 
to the theoretical value of 22, both at the same Mach number 

c v wing chord at span station v 

d v scale factor which has the following values: 

0 . 06 l for v =1 
.234 for v =2 
.381 for v = 3 

.320 for v ' = 4' ' 

Equation ( 2 ) can be written in an alternate form that gives H v 
in terms of wing geometry parameters- that are more significant; thus 


Hy — dy 

where 

K av ratio of average section lift-curve slope to 2 rt/|B, both at the 
same Mach number 

spanwise distribution of section lift— curve slope for a given 
av . Mach number 

C V 

- — spanwise distribution of the wing chord 
c av 

aspect— ratio parameter 

*av 

The term 1/(k/k & v) ( c v/ c av) of equation ( 3 ) gives the effective aero- 
dynamic taper of a wing. The distribution of k/k &v may vary with Mach 
number, particularly at transonic speeds (e.g., due to spanwise variation 
of airfoil section). However, since this distribution contributes to , 
taper effect, the primary change will be in the loading distribution 
rather than in gross load. 

With H v determined from equation (2) or ( 3 ), the values of , 

16 in all, can be found in figure 1 where a^ is given as a function 
of H v for various values of Ap. 

For the case of straight— tapered wings with arbitrary section lift- 
curve— slope distribution for which the chord distribution is specified 
fcy taper ratio, evaluation of equation ( 3 ) is facilitated through use of 


3A 

K av 


( K v/ K av) ( c v/ c av ) 


(3) 
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figure 2, which gives k^/PA for each of the four span stations as a 
function of taper ratio. 


Evaluation of ccy 


The symmetric angle-of-attack distributions most commonly encountered 
are those resulting from .constant angle of attack, twist and camber, and 
the effective angle of attack due to flap deflection. The first two cases 
are the subjects of reference 1. The effective angle-of-attack distribu- 
tions due to flap deflection are the subject of the present report. 

To minimize the computations involved, it is- convenient to consider 
the cases of inboard and outboard flaps independently. The inboard flaps 
are measured from the wing midspan outboard, and the outboard flaps are 
measured from the wing tip inboard. Three flap spans for inboard and 
outboard flaps are defined in appendix A and the equivalent angle-of- 
attack distribution for unit flap deflection ay/S-,, for each case is 
listed in table Ak. The values of ccy/S-j^ given in table Ah are for 
(cp/c) = l. The results can be easily modified for the case of (cp/c) <1 
as will be shown later. The case of full wing-span flaps is easily con- 
sidered since this is the case of constant angle of attack which results 
in additional-^type loading. 


APPLICATION OF METHOD 


For the cases of symmetric angle— of-at tack distributions resulting 
from flap deflection, it is possible to form the set of simultaneous 
equations which are required for the solution of loading distribution for 
an arbitrary plan form. With the loading known, integration formulas can 
be given to determine aerodynamic coefficients. 

The loading distribution G n , determined from the solutions of the 
simultaneous equations, is a function of a vn , which has been shown in a 
preceding section to be a function of wing geometry, compressibility, and 
section lift— curve slope. The aerodynamic coefficients are obtained by 
integrations of the loading distribution and, therefore, will also be 
functions of wing geometry, compressibility, and section lift— curve— slope 
parameters. 

In the following sections, application of the method to obtain the 
general solution for arbitrary spanwise distribution of wing chord is 
outlined and solutions are presented for the special case of wings having 
straight taper. . 
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General Solution 


Simultaneous equations for the loading .- The ay n values to -he 
substituted, in the equations that follow are obtained from figure 1 and 
table I 4 with values of H v given by equation (2) or (3). 

1. Deflected inboard flaps : The flap .spans measured from the 

wing midspan outboard are given by q f . The set of simulta- 
neous equations for symmetric spanwise loading due to deflec- 
tion of any of the three following flap spans can be obtained 
from the appropriate set of the following equations: 


= aixGi+a 12 G2+a 13 G3+a 14 G 4 
= a 2 i5i+a 22 G 2 +a 23 S 3 +a 24 5 4 
= a 31 Gi+ag 2 G 2 +a 33 5 3 +a 34 G 4 

= a 4 lA +a 42 1 % +a 4 3 ^3 +a 4 4 ^4 

where 

G n = G n/ 5 i 

For the case of = 1; a v /6x is unity over the span of the 
wing (same condition as solution for additional loading). 

2. Deflected outboard flaps : The flap spans measured from the 

wing tip inboard are given by qp # The set of simultaneous 
equations for symmetric loading due to deflection of any of 
the three following flap spans can be\ obtained from the appro- 
priate set of the following equations: 



Case 

I 

11 

ill 

3f 

0..195 

0.556 

0.832 

ai 

Si 

-0.017 

0.032 

-o.o 4 i 

a 2 . 

&7 

.029 

-.003 

1.021 

81 

-.014 

.994 

.955 

a 4 

Si 

.988 

.976 

i.o 4 o 


Values of a yn beyond the scope of figure 1 and table I can be obtained 
from equation (B2), which gives the linear asymptotes of the a Vn 
function. 
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Case 

IV 

V 

VI 

V 

0.169 

0.444 

0.805 

Oi 

1.04l, 

0.968 

1.017 

&i 




og 

-.021 

1.003 

.972 

Si 


SSa 

.045 

.006 

1.014 

Si 



2k . 

61 

-.040 

.024 

.012 


- a nGi+a 1 2G2+a 13 G3+a 14 G 4 


— 8-2 iG 1+0,22^2“^ 2 3G3+824 ^4 
= a 3iGi+a3 2 G 2 +a 33 G3+a 34 G 4 


- a 41 G 1 +a 42 G 2 +a 43 G 3 +a 44 G 4 




> 


J 


where 


(5) 


G n = Gn/Si 

For the case of = 1# Oy/Si = 1. 


Loading distribution .— The spanwise loading distributions can he 
determined for various flap configurations as follows: 

1. Full wing— chord flaps [ (c^/c) = 1] : The loading is zero at the 

wing tip and is determined at four other span stations of. 

T| = cos (nrt/8) = 0.924, O.707# 0.383# and 0 by solution of the 
appropriate set of simultaneous equations from the previous sec- 
tion. With equation (A9) and tables A6, A7/ and A8, the loading 
can be interpolated at span stations q = O.981, O.831# 0.558# 
and 0.195 for each of the flap spans considered. With these 
given points and the knowledge that the slope of the loading- 
distribution curve is theoretically infinite at the point of 
angle— of-attack discontinuity or flap spanwise end, the loading 
distribution can be faired. For flap spans other than those 
given# the loading distribution of flaps that are neither inboard 
. nor outboard can be obtained from the concept that loadings are 

additive (discussed later). 


2. Constant fraction of wing— chord flaps [ (c f /c) = constant] : 

The spanwise loading distribution is equal to the product of the 
loading distribution due to full wing-chord flaps and the lift- 
effectiveness parameter da/d&. The factor da/dS is a function 
of the ratio of flap chord to wing chord c-f/c. Values of da/dS 
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are presented in figure 3 , which is reproduced from figure 18 
of reference 4. Although reference 4 from which figure 3 was 
taken limits the Mach number range to values less than 0.2, this 
limitation is believed to be too restrictive since linearized 
theory indicates that da/d 6 is unaffected by compressibility 
for the two-dimensional wing. As indicated in reference 3 , 
da/d 5 is strongly affected by aspect ratio if the aspect ratio 
(or pA) becomes very small; however, comparisons with experiment 
indicate the values of da/d5 from figure 3 appear to be valid 
for J3A > 2. This limit to which da/d5 can be used will be 
considered further in the section Discussion when comparisons 
with experimental results are made. 

3 . Arbitrary spanwise distribution of flap chord [ (c f /c )= variable] : 

The flap can be divided into several parts each- having constant 
da/d 6 and the load distribution due to each part determined. 

The total distribution is then the sum of these individual load 
distributions. 


Lift coefficient .— The lift coefficients due to flap deflection for 
the same flap configurations previously discussed can be found as follows: 

1.* Full wing-chord flaps [(cf/c) = l]: The spanwise loading defined 

at only four points cannot be integrated with sufficient accuracy 
by means of conventional integration formulas. In appendix A, a 
special integration formula is developed that applies to each flap 
span. Equation (a 6 ) can be written as 


PC L6 x 



lllGi + II2G2 + ^ 3^3 + ® 4^) 


( 6 ) 


where, for each of the cases of equations (^) and ( 5 ) ; the hn 
values are given by 



Inboard flaps 

Outboard flaps 

Full- 

span 

Case 

I 

II 

III 

IV 

V 

VI 

flaps 

V 

0,195 

0.556 

0.831 

0.169 

0.444 

0.805 

1.000 

hi 

0.299 

0.299 

0.300 

0.302 

0.301 

0.301 

0.301 

h 2 

.544 

.544 

.555 

.561 

.556 

.556 

.555 


.725 

.725 

.725 

.734 

.728 

.726 

.726 

K 

.392 

.392 

.392 

.405 

.395 

.393 

.393 
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2. Constant fraction of wing— chord flaps [(c^/c) = constant] : 
The flap effectiveness is given by 




to 

dS N K av 


(7) 


3. Arbitrary spanwise distribution of flap chord [(c^/c) = variable] : 

Flaps for which c^/c varies ' spanwise on the wing can be con- 
sidered as an equivalent wing— twist distribution. The effective 
symmetric twist of the wing is given by 



( 8 ) 


where da/dS 
tion (8) is a 
ted by use of 
sidered to be 
loading. 


is now a function of spanwise position. If equa— 
continuous distribution such that it can be plot- 
four semispan points, then the wing can be con- 
twisted and solutions found directly as for basic 


When equation (8) is discontinuous, the angle of attack can 
be divided into spanwise steps of constant angle of attack and 
the total lift can be found by the summation of the lift due to 
each spanwise step. The lift of a spanwise step is obtained from 


a curve of lift coefficient 


as a function of the extent 


of flap span from the wing root outboard. Such a curve can be 
obtained from the solutions of equation (6). As previously 
indicated for = 1 


The lift coefficient due to the twist given by equation (8) 
can be obtained by a method other than the step method from the 
Integral given by 


Taa 


a ( l %x A 


av 


av 


J d6 


cLri 


dq' 


(9) 


which can be integrated numerically by taking the graphical 
slopes of P c L5i A av * 
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Spanwise center of pressure, induced drag, and pitching moment .— 

The integration formulas of loading due to flaps required to obtain span- 
wise center of pressure and induced drag are not developed here. However, 
the integration formulas given in reference 1 to obtain these characteris- 
tics can be used with acceptable accuracy for the case of loading due to 
flap deflection. 

No rigorous procedure for determining pitching moments can be devel- 
oped from this method. As noted in reference 1, however, a good approxi- 
mation for a wing without flaps can be obtained if the wing sweep is large 
since moments due to shifts of spanwise load overshadow those due to shifts 
of chordwise load. Where the effect of flaps is to be considered, as in 
the present case, it appears unlikely that the shifts in chordwise loads 
can be safely neglected. Therefore an estimation of pitching moments due 
to flaps should not only include the moment due to spanwise load redistri- 
bution, which this method will give accurately from the longitudinal 
moment of the center of load, but also should include an estimation of the 
moment due to chordwise redistribution. 


Additional considerations,— Several items pertaining to the usage 
of the method, should be considered. 


1. Additive nature of loading and spanwise angle of attack: 

The linear relation between angle of attack and loading distri- 
bution of equation (l) states that all loadings are additive if 
the respective angle-of-attack distributions are added. 

Thus , the aerodynamic coefficients that result from the 
integration of .the loading distribution are additive if they 
depend linearly on the loading distribution. Hence, lift coeffi- 
cient, rolling-moment coefficient, and pitching moment are addi- 
tive; whereas induced drag, spanwise center of pressure, and 
aerodynamic center are not. 


The additive concept is very useful for the determination 
of loading due to flap. deflection. The loading due to an arbi- 
trary span flap having arbitrary position on the wing . semi span 
is simply found by adding or subtracting known loadings due to 
flaps at other locations on the wing. For example, if the span- 
wise ends of the flaps are at the same span stations, then for 
full-chord flaps 


G inboard “ G tu=i G outboard 


(%Xw. rf =0L a-( C O 


inboard 


outboard 


= 1 - TJ f 

inboard outboard 



( 10 ) 
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2. Differential flap angles : The effect of a differential between 

flap angles can be taken into account by considering C-^ of 
each wing panel as one— half the results of the symmetric ^case. 

The total wing lift is then the sum. of the products of Cj^/2 

and the angle of deflection of each flap. Although the total 
wing lift can be found by this procedure, the loading distri- 
bution cannot be found exactly by the products of the symmet- 
ric g/ 5 and the deflection of each flap. However, the load- 
ing distribution will be sufficiently accurate for most purposes 
since this procedure neglects only the small change due to the 
induced effects of the difference in flap deflection on each 
opposite wing panel. 

3. Flap angles measured perpendicular to the hinge line : The rela- 

tionship between a flap angle measured in a plane perpendicular 
to the flap hinge line and that measured in a plane parallel to 
the plane of symmetry is given by 

tan 5 = — 6 (11) 

cos Aj. 

where is the sweep angle of the flap hinge line, and 5 

is the angle measured perpendicular to the hinge line. Equa- 
tion (ll) can he inserted into any of the equations involving 
5 to put the equations in terms of 6. 

For constant fraction of wing— chord flaps on straight— tapered 
wings, A f is given by 


tan Aj, = tan A— ^A) ~ ( c fA)] / i_ x \ 

a vi+xy 


Solutions for Straight— Tapered Wings 

Charts of aerodynamic characteristics for straight— tapered wings can 
he presented in terms of parameters involving geometry, compressibility. 
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and section lift-curve slope. 5 These charts provide a ready means of 
obtaining data directly. . 

Equation (.4) has been vised to evaluate the loading distribution 
G/8j. = (i/Sjl) (cjc/2b) for a range of wing parameters and flap spans, 
and the results are plotted in figure 4. The loading presented in fig- 
ure 4 gives the values due to full wing— chord flaps. For results of 
other flap configurations, the procedure given in the loading distribu- 
tion section of the general solution can be used. 

With the loading distribution determined from equation (4), the appli- 
cation of equation 
in the form |3C T 

L &i 

are plotted against extent of flap span from the wing root outboard for a 
range of wing parameters. 

As presented, figure 5 gives directly the lift due to full wing— chord 
inboard flaps for flap spans measured from the plane of symmetry outboard. 
The lift due to outboard flaps for flap spans measured from the wing tip 
inboard can be obtained from figure' 5 hy use of the relations of equa- 
tion (10). For full wing-chord flaps located arbitrarily on the wing 
semispan, the lift can be obtained from figure 5 as indicated in the folr- 



(6; provided lift coefficients due to flap deflection 
/k &v for several flap spans. In figure 5, these values 


throughout the figures, k &v is the constant spanwise— section lift-curve 
slope or the average of a small variation. For large spanwise varia- 
tions of k that follow the function given in equation (B5) developed 
in appendix B, the parameters PA/k and X can be replaced by the 

parameters j and ^ X,- respectively. 

For large spanwise variations of k that do not follow the curve of 
equation (B5), the simultaneous equations for the general solution can 
be solved for arbitrary distributions of k. Values of H v can be 
obtained conveniently from figure 2. 
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With the full wing— chord ■values given above, the lift due to constant 
fraction of wing— chord flaps and flaps of arbitrary spanwise chord distri- 
bution can be found through use of equations (7) and (9) with the da/dS 
values of figure 3» 


APPLICATION OF THE a Vn COEFFICIENTS TO THE PREDICTION 
OF DOWNWASH FROM AKBITRAEY WINGS 


In the subject report and in reference 1, it is shown how various 
forms of symmetric loading can be found quickly through use of simplified 
lining— 1 surface theory and values of the a vrL coefficients which are pre- 
sented in graphical form. It can also be shown that if the symmetric load 
distribution is known, from any computational method or from experiment, 
then the same ay-^ coefficients can be used to find corresponding values 
of the downwash in the center of the downwash wake. It is the purpose of 
this section of the report to make available the method for computing 
downwash without attempting to explore the possibilities or limitations 
of the method. In view of the importance of being able to predict down- 
wash, it is considered that inclusion of this material herein is justifi- 
able although not directly related to the main purpose of the report. It 
is important to remember that the method can be used for any, case of sym- 
metric loading and is not restricted to that discussed in this report. 

As noted previously, the a y n coefficients provide a very simple 
means of finding the downwash in the center of the downwash wake of the 
wing due to any symmetric spanwise loading distribution for any sweep 
a ngls. The downwash can be found at specific points in the region between 
the wing tips from near the quarter-chord line to an infinite distance 
downstream. 

Equation (l) was derived (reference 1) to find the downwash at the 
effective three— quarter— chord line for which the distance downstream from 
the quarter-chord line is k(c/2). The general expression for the distance 
downstream from the quarter— chord line is 


x — (3y tan Ap 

where x is the longitudinal distance measured in the horizontal p i a me 
from the wing-root quarter-chord point to an arbitrary downwash point. 
Substituting this general expression in equation (2) gives 


H v = 


U/P) - T| v tan Ap 


(13) 
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Equation (l) (for ■values of the parameters i/$, r| v , and tan A^) gives 
the downwash angle as P 



(l4) 


where the a Vn coefficients are obtained from figure 1 with the H v 
values given by equation (13). 

It should be noticed that, because the chordwise loading is concen- 
trated at the one-quarter— chord line, equation (13) is independent of 
aspect ratio and spanwise distribution of wing chord, while from equa- 
tion (l4) it can be seen that downwash depends on these wing geometric 
values only through the loading distribution G n . 

Equation (l4) will give the downwash angle due to arbitrary symmet- 
ric loading distribution for all values of > q v tan Ap. The downwash 
can be found for the specific span stations corresponding to V= 1, 2, 3> 4, 
or q v = cos ~ = 0.924, 0.,707, O.383, and 0, respectively. If compressi- 
bility effects are included, the loading distribution must be that at the 
Mach number considered. 

For convenience, equation (l4) is rewritten such that downwash angle 
and loading distribution are in terms of unit lift coefficient. Then 


k 



where c x c / c L c av is " thie s P an loading distribution for unit lift. Both 

a „ and c-,c/c T c 0 „ are affected by compressibility, being functions of 
vn. Jj civ 

the parameters 0A/K av and Ap. 

Equation (l4) or (15) gives the downwash, in the vertical center of 
the wake, due to the integrated effect of the loading of a swept vortex 
and the trailing vortex sheet extending backwards to infinity. ' Any roll- 
ing of the trailing vortex sheet into two trails is neglected. Since 
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either of the equations are for arbitrary loading distribution, the 
loading coefficient can be set equal to a constant to duplicate uniform 
loading and the downwash results obtained for the case of the trailing 
sheet rolled into two trails. 

With C 2 c/c L c aT = 1, equation (15) becomes 


pc lA ar' 


where 

4’ 


TJ* 

and b' (as used in the past from line— theory concepts) is the span 
between the two trails given by the integrated loading divided by the 
loading at the middle of the wing, or 


{c l c/C -L c ^\ =0 


DISCUSSION 



( PA 

Ik. 

\ av 


k 

L' 

n=i 


Vn 


(16) 


- 


b»/2 

y_ 

b*/2 


Influence of Wing Geometry on Flap Effectiveness 
for Straight— Tapered Wings 


To show the effect of wing geometry and compressibility on lift due 
to inboard flap deflection, the data of figure 5 have been cross-plotted 
and presented in figure 6. These effects are similar to those on wing 
lift-curve slope shown in reference 1. Figure 6(a) compares the present 
theory with that of reference 3 and shows, as might be expected, that 
both theories give essentially the same results for -M- < 1. Figure 6(b) 

av 

shows the influence of taper ratio on flap effectiveness for swept wings. 
Figure 6(c) gives the relative effects of the wing parameters of sweep, 
aspect ratio, and taper ratio. 
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Figure 7 shows comparative effects of inboard and outboard flaps on 
swept wings and shows the loss of flap effectiveness due to sweepback is 
less for inboard flaps. The relative effects indicated by figures 6 
and 7 apply equally well for constant fraction, of wing— chord flaps since 
the lift is only reduced by a constant spanwise factor da/d8. 

Limitations of the da/dS Concept and 
Compressibility Considerations 


A minimum value of aspect ratio for which two-dimensional values of 
da/dS can be used with accuracy has been determined for use with the 
subject report from principles of lifting— surface theory and from experi- 
ment. An analysis based on lifting-6urface theory indicates that the 
da/d6 of a wing does not vary appreciably with aspect ratio until aspect 
ratio becomes less than approximately 2} for A < 2, da/d6 increases 
rapidly with decreasing aspect ratio until at A = 0, da/d6 = 1'. The 
discussion of reference 2 concludes that only a lifting-surface theory 
that accounts for changes in chord loading due to flap deflection will 
give accurate loading predictions for low-aspect— ratio wings. 

Low-speed experimental values of da/d&, from a 7— by 10-foot wind 
tunnel (reference 5), are plotted in figure 8 as a function of aspect 
ratio for an unswept and a swept— back wing. Figure 8 indicates the 
order of variation of da/d& with aspect ratio.' 

The compressibility considerations in the present report are the 
same as those of reference 1, 2, or 3i namely, that an incompressible 
theory can be made into a compressible theory at subsonic speeds subject 
to the limitations of the three-^dimensional, linearized, compressible- 
flow equation. As discussed in these references, compressibility can be 
taken into account by substituting for the geometric values of A, A, 
and b, the parameters PA, Ap, and 0b, respectively. 

From two-dimensional theory it can be shown that the parameter 
da/d6 is independent of Mach number . In the present report da/d8 is 
also considered independent of Mach number and, further, is considered 
independent of aspect ratio for A > 2, as discussed earlier. However, 
from the linearized compressible— flow theory the aspect ratio becomes an 
effective aspect ratio given by 0A which approaches zero as M 1.0. 
Thus the use of two-dimensional da/dS is restricted to the Mach number 
range for which 0A > 2. 

Figure 9 shows the flap-effectiveness comparison, between theory and 
experiment for a subsonic range of speed for two wings. The theoretical 
values were obtained using da/d& from figure 3 and assuming it inde- 
pendent of Mach number. • The curves are terminated at the Mach n umb er 
where 0A = 2. The parallelism of the theoretical and experimental curves 
of figure 9 indicates that d<x/d& can be used with the assumption that 
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it remains independent of Mach number to values of PA of the order of 
two. 


In view of the foregoing discussion, it is believed that two- 
dimensional values of da/d 5 can be considered independent of M as 
long as PA ^ 2 in the application of the subject procedure. 


Comparison of Theory with Experiment 


Experimental values of flap effectiveness are compared with theoreti- 
cally predicted values in a correlation diagram given by figure 10. 

Included are the results of a wide range of plan forms. Sketches of the, 
plan forms and flaps are drawn about the points of correlation. The theory 
makes use of the curve of figure 3 giving da/d5 for a sealed-gap flap 
over a range of deflection of ±10°. As seen in figure 10, the theory is 
in fair agreement with experiment. 

In general, the theoretical' values of figure 10 are too low and 
center approximately about the minus 10— percent line of the figure. A simi- 
lar effect is shown in figures 8 and 9. This almost constant difference 
between the theoretical and experimental results is possibly due to the 
da/dS values that were used in obtaining the theoretical results and 
which were obtained from the —10° to 10° sealed— gap curve of figure 3« 

The theoretical flap effectiveness of the low-aspect— ratio wings 
having partial-wing-chord flaps, given in figure 10, fall below the 
minus 10-percent line. Part of this discrepancy is presumably due to 
using two-dimensional values of da/d6 since these wings approached the 
lower limit of aspect ratio for which the da/d& concept has been judged 
acceptable. 

Comparison of experimental results with theory for a Mach number 
range is shown in figure 9. This figure shows that compressibility effects 
are accurately predicted up to Mach numbers corresponding to values of 
3A -> 2. 


CONCLUDING EEMARKS 


The determination of symmetric loading due to flap deflection for 
arbitrary wings is shown to be obtained by the solution of four simulta- 
neous equations. The coefficients of the simultaneous equations are pre- 
sented in charts as functions of parameters that include wing geometry, 
compressibility, and section lift— curve slope. Integration formulae are 
given for the loading distribution. With these the flap effectiveness 
can be found. 1 
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For the special case of straight— tapered wings, the loading distribu- 
tions and values of flap effectiveness are given in chart form for a range 
of wing plan forms. These charts indicate that the effects of plan form 
on the loading distribution and lift due to flap deflection are similar 
to that on additional loading. 

On the basis of experimental comparisons, it is concluded that the 
method of the subject report can predict the flap effectiveness of wings 
at subsonic speeds with good accuracy. 


Ames Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Sept. 22, 1950. 
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APPEND IX' A ' ' 

DEVELOPMENT OF EQUATIONS FOR DETERMINATION OF SYMMETRIC 
SPAN LOADING AND LIFT DUE TO FLAP DEFLECTION 


Altered Boundary Conditions for Flap Deflection 

In the following development the flap chord is considered equal in 
length to the wing chord; thus a deflected flap results in a discontin- 
uity in the spanwise distribution of angle of attack. Spanwise loading 
for a discontinuous distribution of angle of attack can be found through 
the simultaneous solution method of equation (l). The number of equa- 
tions, however, must be greater than four to obtain accurate results and 
then the computation efforts become excessive. (Reference 1 gives the 
method for an arbitrary number of equations.) However, an alternate pro- 
cedure, originally developed in reference 2, provides sufficient accuracy 
with a reasonable number of equations and this procedure is the basis of 
the following development. 

The zero-aspect— ratio theory of reference 3 shows that for a given 
sweep angle, as aspect ratio approaches zero, spanwise loading becomes 
independent of plan form. Thus, the altered boundary conditions (see 
Development of Method) required for the present theory to give spanwise 
loading for given flap spans can be determined from the loading given 
by reference 3> which applies to all plan forms. These altered boundary 
conditions are obtained from equation (l) using the a • values for zero 

aspect ratio and the loading distribution for zero-aspect— ratio wings 
given by reference 3 for given flap spans.. 

The zero-aspect— ratio values of a^ are given by figure 1 for 

Hy = 0 or are given by the B Vn coefficients of reference 1. These 

a., coefficients are tabulated as follows: 

vn 


TABLE A1 

a vn (for m = 7 and A = 0) 


B 

l 

2 

3 

4 

.1 

10.4524 

-2.0720 

0 

-0.2242 

2 

-3.8284 

5.6568 

-1.8284 

0 

3 

0 

-2 . 3888 

4.3296 

-3.1548 

4 

-.2928 

0 

-1.7072 

4.0000 
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With these a vn coefficients and equation (l), symmetric spanwise 

loading can he found for zero-aspect— ratio wings. As a comment on the 
accuracy of the simplified lifting-surface theory for determining addi- 
tional loading with m = 7, the solution of equation (l) with A = 0, 
a vn v& l ues gave, to the accuracy of three decimal places, the elliptic 
loading distribution characteristic of zero-aspect— ratio wings. 

The flap end is arbitrarily chosen for the present theory as the 
mean value of the spanwise trigonometric coordinates of the points at 
which the boundary condition is applied and which bracket the discon- 
tinuity or flap spanwise end. Thus for m = 7> three flap spans can 
be defined for both inboard and outboard flaps. Let T)j> be the flap 
span, and 0 the spanwise point at the end of the flap, then 


■p ( £ = cos 0 for inboard flaps 
Tlj. = 1-cos 0 for outboard flaps 
These flap spans are given in the following table: 


TABLE A2 

DEFINED FLAP SPANS 



Inboard 

Outboard 

Case 

I 

II 

III 

IV 

V 

VI 

0 

2* 

5* 

21 

3* 

21 

11 . 

16 

16 

16 

16 

16 

16 ' 

V 

0.1951 

0.5556 

0.8315 

0.1685 

0.4444 

0.8049 


For the flap spans listed in table A2, the exact span loading 

distribution for zero-aspect ratio can be found from reference 3« 

With the a values listed in table A1 and the exact values of 
vn 

G x , G 2 , G 3 , and G 4 determined from reference 3> equation (l) gives 
the boundary condition or twist required for the present method to 
produce the loading distribution for each case listed in table A2, 
that is. 
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2a = 10.4524 2i -3.8284 2§- -0.2928 2± 

5i &x &i 8 X 


2s. = -2.0720 2i +5.6568 2i. -2.3888 — 

81 81 8 X 5 x 


2a = _i.82.84 2^ +4.3296 -1.7022 2d 


2t = -0.2242 2i -3.1548 2a +4.0000 2t 
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The spanwise loading distribution due to flap deflection is given 
by (reference 3): 


Inboard, flaps 


Siri. = S i n cp _ I 

81 rt 


29 sin cp + (cos cp — cos 0) Zn 


(cos 9 + cos 0) Zn 


cos 


0+9 


cos 




Outboard flaps 
~ G(q>) 

- i8x 


outboard 


■ sln * - [Itj 


0+9 

— n — 2- 


sin 


0-9 


inboard 




(A2) 


With equation (A2), the spanwise loading Gx, G&, G 3 , and G 4 , at span 

K 


it rt 3* 

“ 8' 4' 8 J 

tan be tabulated for each of the cases given in table A2 as follows: 


stations 9 = *, and - or q = 0.9239, 0.7071, O.3827, and 0 

8 4 8 2 
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TABLE A 3 

SPANWISE LOADING 'FOR A = 0 


81 

Case I 

.11 

III 

IV 

V 

VI 

Gi 

57 

0.0505 

0.1524 

O.2586 

0.1241 

0.2303 

0.3322 

G2 

57 

.1106 

.3467 

.6432 

.0639 

.3604 

.5965 

g 3 

Si 

.2062 

.6903 

. 8714 

.0525 

.2336 

.7177 

g 4 

Si 

.4125 

.7970 

.9618 

.0382 

.2030 

.5875 


The twist distribution required for each case is obtained with equa- 
tion (Al) and table A 3 , and is summarized as follows: • 


TABLE A 4 


TWIST DISTRIBUTION, — 

81 



Inboard 

Outboard 

Case 

I 

II 

III 

IV 

V 

VI 

e 

ll 

2l 

21 

21 

5 * 

11 


16 

16 

16 

< 16 

16 

16 

n f 

0. 1951 

0.5556 

0.8315 

0.1685 

0.4444 

0.8049 

<Xi 

Si 

0170 

.0318 

04i4 

1.04i4 

.9682 

1.0170 

a 2 

57 

.0285 

-.0034 

1.0210 

-.0210 

1.0034 

.9715 

\m 

-.0136 

.9941 

.9548 

.0452 

.0059 

1.0136 

a* 

57 

.9884 

.9763 

l.o4o4 

-. o4o4 

.0237 

.0116 


With the twist distribution given by table A 4 , equation (l) can be 
used to solve for spanwise loading due to flaps for any of the six cases 
and any aspect ratio. 
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Lift Due to Flaps 


The loading distribution resulting from flap deflection is distorted 
to such an extent that the usual quadrature formulas are not sufficiently 
accurate to integrate the loading for m = 7. However, the following devel- 
opment results in a method which will give the lift due to flaps as a func- 
tion of the four lift coefficients: 

Lift coefficient is given by the integral 




G(n)dn 


A 



G(9) sin cpdcp 


If it is assumed the loading series 

00 

G(q>) = a 4i sin Hi 9 

then 




(A3) 

If the loading series is expanded for 9 = > an<i gt an ^ 

Gx,. G2 , G 3 , and G 4 are obtained in a series of a^'s, then multiplying 
each series through by the first coefficient of the series and summing the 
G*s gives 

a l = h (0.3827 Gx + 0.7071 G2 + 0.9239 0 3 + vp) + a i5 -a 17 +a 3i- a 33 + • • • 

, (a4) 


The high harmonic coefficients can be put as factors of the G n . 
With equation (Ah), the lift coefficient of equation (A3) becomes 


Cl = A 


0. 3827rt 
% 




1 + ^-^ 17 +a 31 


1.5068 G 




Gx + 


0.70717T 


c 




1.5068 G 2 


) 


+ 


0.9239^ 




1.5068 g 2 


1 + a lFS~ a 17+ a 31~ a 33 ^ G + £ f 1 + a lS~ a 17 +a 31~ a 33 


¥V2 


1.5068 G„ 


) 


G. 


(A5) 
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If h n is defined as the coefficients of G n , then the lift coef- 
ficient is given by 


Cl - A (hiGi + h 2 G 2 + h 3 G 3 + h 4 G 4 ) 


(A C) 


Reference 3 also gives loading due to flaps in series form from which 
the ratio of (a 15 — > a i 7 +a 3 i -a 33 ) G n can eva l ua 'ted. It is expected 

that this ratio will not vary appreciably with aspect ratio. These high 
harmonic coefficients are small but are not negligible for loading due to 
flaps. The values from reference 3 are given by the following equa- 

tions (for odd values of a^): 


For inboard flaps 


Pi = 1 


ai = (n-2 9 + sin 29) 

n 


pi > 1 


hb x 


a^i = ]tp 1 (^ 1 g_i) sin 0 cos ^1 0 “ cos 0 sin ^l 0 ) 

For outboard flaps (for the same 9) 


(A7) 


Pi = 1 


( a i) ^ • = &i - (ai) 

outboard inboard 


Pi > 1 


(a+ij = - (aji x ) 

outboard inboard 


With equations (A5) and (A7) and the values of table A3, the 
values for each, case can be tabulated: 
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TABLE A 5 

LOADING INTEGRATION FACTORS, 


Case 

I 

II 

III 

IV 

V 

VI 

- hi 

0.2991 

0.2994 

O.2999 

0.3020 

0.3014 

0.3008 


. 55^1 

. .5544 

.5549 

.5606 

.5563 

.5556 

h 3 

.7248 

.7250 

.7252 

.7339 

.7275 

.7259 

h 4 

.3922 

.3921 

.3922 

.4050 

.3950 

.3930 


Interpolation Function for Spanwise Loading Distribution 


The method of this report gives spanwise loading due to flap deflec- 
tion at four span stations. Since the complete loading distribution is 
not sufficiently defined by these values, the values of loading at other 
span stations cannot be interpolated accurately by direct use of the inter- 
polation equation and table of reference 1 . 

The direct use of the interpolation table of reference 1 can be 
obtained by converting the present loading due to flap deflection to a 
loading approximating additional loading distribution. The conversion . 
factors are obtained from values of the zero-aspect— ratio theory and 
given by the ratio of equation (A 2 ) to sin q> n for each of the four given 

span stations. Define 


r = ( G / Si )a= 0 
sin cp 


(A 8 ) 


then, with table A 3 and equation (A 8 ), values of R n have been computed 
and are tabulated for each of the flap span cases. 
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TABLE A 6 



Inboard 

Outboard 

Case 

I 

II 

III 

IV 

V 

VI 

n 


0.195 

0.556 

O.832 

0.169 

0.444 

0.805 

1 

0.924 

0.1318 

0.3981 

0.6757 

0.3243 

0.6019 

0.8682 

2 

.707 

.1564 

.4903 

.9096 

.9041 

.5097 

.8436 


.383 

.2232 

.7471 

.9431 

.5685 

.2529 

.7768 

4 

0 

.4125 

• 7970 

.. 9618 

.0382 

.2030 

.5875 


where, for qp = 1, B n = 1. 


The interpolation table of reference 1 gives 


TABLE A 7 

INTERPOLATION FACTORS e^ 

I . . - 1. * A I v\. c I. /%. C-C.C /l .1 A 1 AC.1 

ERRATA NO. 1 
NACA TN 2278 


THEORETICAL SYMMETRIC SPAN LOADING DUE TO FLAP DEFLECTION 
FOR WINGS OF ARBITRARY PLAN FORM AT SUBSONIC SPEEDS 
By John DeYoung 
January 1951 


The following equation should be inserted at the bottom of page 29? 


4 



n=i 


NACA-Langley - 3-18-53 - 1000 
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where equations (A2) and (A8) provide values of R k which are tabulated 
as follows: 


TABLE A8 
*k 


•k 

"'^v*Case 

I 

II 

III 

IV 

V 

VI 

1 

2 

0.981 

0.1246 

0.3826 

0.6090 

0.3911 

0.6176 

0. 875^ 

3 

2 

OJ 

00 

CO 

• 

. 1421 

.4304 

.8090 

.1910 

.5695 

.8580 

2 

.556 

.1799 

.6253 

.9314 

.0685 

.3746 

.8201 

2 

2 

.195 

.3324 

.7863 

.9556 

.0442 

.2136 

.6675 


With the R k values of table A8 and the loading ratios determined' 

R k 

from equation (A9) , the loading G k can be found at span stations 
k = ^ , and ^ (q = 0.981, O.832, 0.55 6> and 0.195> respectively.). 


\ 
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APPENDIX B 
AIDS TO THE METHOD 

Linear Asymptote of a 

vn 

For large values of Hy , the ay n functions become linearly propor- 
tional to Hy (see, e.g., fig. 1). Since this linear characteristic appears 
at relatively low values of the simple linear relation between ay n 

and Hy is quite usable. 

The L v ^* function of appendix A of reference 1 is multiplied by 
b/c v and the product is linearized with respect to b/c v . The result 
gives 
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With the values of equation (Bl) substituted into equation (A37) of 

reference 1, the values of a vn _ for arbitrary sweep angle and for high 

values of H v are obtained. For m = 7# the following equation (B2) 

gives values for a.^ as 

vn 


= ( + 1.452 tan AJ Hi + 7.915 - 1.751 sin A - 

\cos A / 

0.003 tan £ + 0.088 + _ 

2 \ tan A J 


0.052 (l 1 ^ °- 01 77 W A-^V 00l4 

^ tan A ' 


■^)+ o.oi4 ^ ^ 1+ °- 


1716 tan 2 A - l' 
tan A - 


a 12 =( Q- 8 ^ 2 - 3.793 tan A JH X - 3.100 + 1.6l0 sin A + 
\ cos A / 


0.008 tan ^ -0.048 


1 + 0.0016 tan 2 A - 1 
tan A 


0.034 


1 + 0.0177 tan 2 A - 1 
tan A 


- 0.052 


1 + 0.1716 tan 2 A - 1 
tan A 


( ~° ,g ^ 6 + 2.476 tan a') H x - 0.014 - 1.727 sin A - 
\ cos A / 


0.032 tan ^ + 0.037 ( — - 0,0016 - 

d K tan A 


=0 


. 1251 ^ 


tan A 


1 + 0.1716 tan 2 A _ 1 


tan A 



t 
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0.084 


cos A 


- — 1.110 tan A ) H x - 0.043 + O.637 sin A - 


0.060 tan ^ - 0.017 ( -+ ■ PP 1 ^ , ,, A 

2 ' tan A ' 


0.048 


/ /rr 


0.0177 tan 2 A - 
tan A 


i) ♦ 0 omf' -l 1 + 0 - 1716 w a ^.1 ) 

' 7 tan A ' 


a 2 i = ( 0,231 + 0.408 tan a) I 2 - 2.140 - 0.223 sin A - 
\ cos A ' 


».004 tan ± - 0.063 ( ^ A - 1 ) + 


tan A 


0.116 6 V 1 .. * 0.0294 + 0 018 / yr 

' tan A ^ ' 


+ 0.0886 tan 2 A — 


tan A 


a 22 = ( °- 9 ^ 1 + 0.267 tan Aj H 2 + 4.640 - 0.491 sin A + 

\ cos A / 

0.011 tan - + O.163 ( ^ 1 - °: . 01 T - 1 ta . ng A ~ A _ 

2 ' tan A • / 

0.063 ( Jl * Wa - * ) 0.068 (A * 0-0886. tan* A- 1 ) 

' tan A 7 ' A ' 


tan A 



3 ^ 


a S3 


( 0 . .2J 5. _ 0.986 tan A) H 2 - 2.226 + 0.821 sin A - 
\cos A J 
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0.042 tan * - 0.107 + °.°m to* A- 1 

2 V tan A 


o.oi+8 ('£- 1 - Q- 0 - 2 .^ - 0.063 Q- . 

\ tan A / \ 


0886 tan 2 A — 


tan A 
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_ 0.06 . 7 . + 0.377 ^ A h 2 + 0.188 - 0.551 sin A - 
cos A / 


0.079 tan | ♦ oW: ) - 
2 V tan A ' / 

0 022 ( ^1 + 0.0294 tan 2 A - lA + 0 . n6 / V 1 + 0.0886 tan^A - l N 
• \ tan A J * \ tan A / 


>(B2) 


*31. 


= ( _ 0,041 _ 0.068 tanA J H a + 0.197 + 0.207 sin A - 
\ cos A / 


0.007 tan 


A _ 0 n6 A/l + 0.1716 tan 2 


A - 1 


V 


tan A 


125 f */l + 0.0886 tan 2 A - l \ + Q 21 ^ ( JjL + 0.1993 tan 2 A - l \ 
\ tan A ' ^ tan A ' 


a = [ 0020 + o.251 tan A) H„ - I.85O - 0.132 sin A + 
32 V cos A J 

*/l + 0.1716 tan 2 A - l' 


' 0.020 tan ^ + 0.302 


tan A 


0 


.082 ( A - 1 . Q « 0886. tan 2 A -7.1.^ _ o.ll6 A . t.Q i. 
\ tan A / \ 


1993 tan A - 1 


tan A 
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0-543 


a M = tan a m 3 

V cos A J , 


+ 0.232 tan A] H 3 + 4.272 - O.388 sin A - 


0.078 tan — — 0.197 \ + °' 1716 ^ A - ^ + 

2 \ tan A J 


0.302 


'l + 0.0886 tan 2 A - 
tan A 


a) + 0.088 


1993 tan 2 A - 
tan A 


= ( OJL56 _ Q 1^28 tan A^ — 1.666 + 0.489 sin A - 
Vcos A J ^ 


0.14? tan — + 0.088 — ± W A ~ 1 ')- 

2 V tan A / 


o.ii g( £ ± ° ’ 0886 tan2 ^.-.- 1 -^ - O.i 
V tan A / 


r l + 0.1993 tan 2 A 
tan A 


a. . = 0.032. ( 1 . + tan A^ BL - 0.224 
4 - L \cos A y * 


a 42 = -0.097 ( — 1 — + tan A ) H4 + 0.354 ( 1 + sin A 
Vcos Ay V 


a, „ = 0.371 (■ — + tan A) H 4 - 3.155 

\cos A / 

a,, = O.696 { — i — + tan A j H 4 + 3-313 — O.687 sin A 

** . V cos A / 



36 


NA.CA TN 2278 


Linear Spanwise Distribution of (Kc) v /(itc) av 

Equation ( 3 ) can be manipulated into a function of various parameters 
involving distributions of wing chord and section lift— curve slope. One 
procedure is as follows: With the condition that the product of section 

lift-curve slope and wing chord varies linearly spanwise, then 


(KC) = 


2b 






(B3) 


and equation ( 3 ) becomes 

• 1 + Kj.x/k 

HV “ dv<M " ) 2tl - -| V (1 - K t X/K r )] 

where A K is the aspect ratio based, on a wing chord equal to (kc). Equa- 
tion (B3) is reduced to terms of two parameters, A K and ( K t/ K r)^ # 
'Expressions of Ak in terms of aspect ratio for straight— tapered wings, 
'and the distribution of section lift— curve slope can be found. 

For straight— tapered wings 


c r (l + X) 

and since (kc) is linear 


Ak - 


2b 

K r c r [! + (K t /K r )X] 


then 


A(l+X) 
K r + K t* 


and equation (B3) becomes 


Pa( i+x) 1 1 + ( K t A r )>. 

K -tx) J j 


(B4) 


Hy - dy 
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The distribution of k for straight— tapered wings is given by 


(kc) 


K\i — 


V = K, 


1 _ H v [l ~ (* t /* r )M - 

1 - ti v ( 1 - X) 


(B5) 


ftA ( 1+X ) 


and 


Equation (BU) is in terms of two parameters given by 

L (K r + K t \) J 

(K^./K r )X. Solutions for' spanwise loading in terms of these two parame- 
ters and are valid for the distributions of section lift— curve slope 

P 

given by equation (B5). Equation (B5) indicates that at X = 1, K v is a 
linear function and at X = 0, k v is a constant. For values of X 
between 0 and 1, k v is a curve in the region between the linear function 
and a constant. 


Equation (b 4) is represented in graphical form by figure 2 but with 

IL. 

the ordinate given by the parameter — ; 

*, (K t A r ),. 


and the abscissa 


In summar y, for straighi^-tapered wings having a linear distribution 

(kc) , the loading and associated aerodynamic characteristics can be pre— 
v 9 

sented for a range of the parameters V and - (K t/ K r ) x - 
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Figure /.- Influence coefficients, a, p , for symmetric spanwise loading plotted as a 
function of the wing-chord-distribution parameter, H„ for values of the sweep 
parameter, A 0 , degrees. 
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Figure l- Continued. 
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. (j) * • 3, n = 2, d 3 m 38i . 
Figure ! .-Continued. 
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Figure /- Continued. 
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(!) y = 3, n* 4, d 3 = .38/ 


Figure l- Continued. 
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Figure l- Continued. 
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Figure /.-Continued. 
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Figure 2r Variation of the geometric parameter, 


Ky Hy 

w 


with taper ratio, A, for 


straight -tapered wings. 



Figure 3.- Variation of lift -effectiveness parameter with flap -chord ratio, 

Cf 

— . Average trailing -edge angle about 10° , M -0.2 . (Curves from 
reference 4.) 
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(a) k^- s 0, Inboard and outboard flaps for all plan forms 
with the trailing edge at or behind the upstream 
line of maximum span (reference 3). 


Figure 4. -Spanwise load distribution, per radian, due to 

symmetric flap deflection for straight - tapered wings. ££. =/ 

c 




.8 .9 


Fraction of wing semi span, 7 
(b) &A = 2 q. inboard flaps. 

\ 

Figure 4. -Continued. 













Fraction of wing semispan ) y 
(d) 10.0) Inboard flaps. 

Figure 4. -Concluded. 
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./ .2 .3 .4 .5 .6 7 .8 .9 10 

Fraction of wing semispan, y 

(o) X - O . 

Figure 5.- Flap lift parameter, — — per radian, for extent of flap span 


from the wing root outboard , 


* /. 



Fraction of wing semispan , rj 
(b) X -- 0.5 . 

Figure 5.- Continued . 
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Fraction of wing semispan, 7 
(c) X s t.O . 

Figure 5. - Concluded . 


(a) Effect of aspect ratio for X - 0.5 . 



Figure 6. - Influence of wing geometry on flap effectiveness (%-=/) 






per radian 
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A 0 , degrees 

p 

(c) Effect of sweep for - 0.4 . 


Figure 6 - Concluded . 
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la) Un swept wing, ‘ O’. 



(b) Swept wing . 

Figure 8. - Variation of experimentally determined da/ds with 
aspect ratio for wings having flaps of 30% of the wing 
chord normal to the c/4 tine. - t.O, x = 0.5 . 





Mach number, 


MACA TN 2 



Reynolds number = 2x10 s at M*0.2; 4.8x10 s at M s 0.6 
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O .5 10 15 2.0 2.5 3.0 

Experimental C L s t per radian 


Figure iO -Correlation of theoretical and low- speed experimental f/ap- 
effectiveness due to symmetrically deflected flaps for various 
plan forms: Reynolds number from / to 3xKf. M<0.2. 


NACA -Langley - 1-16-51 - 1000 




